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j§ : ABSTRACT 

> : 

^ ■ We present a non-perturbative solution of large iV matrix models modified by terms 

of the form g(Tr $ 4 ) 2 , which add microscopic wormholes to the random surface geometry. 
For g < gt the sum over surfaces is in the same universality class as the g = theory, and 
the string susceptibility exponent is reproduced by the conventional Liouville interaction 
~ e a+< ^. For g = gt we find a different universality class, and the string susceptibility 
exponent agrees for any genus with Liouville theory where the interaction term is dressed 
by the other branch, e a ~^. This allows us to define a double-scaling limit of the g = gt 
theory. We also consider matrix models modified by terms of the form gO 2 , where O is 
a scaling operator. A fine-tuning of g produces a change in this operator's gravitational 
dimension which is, again, in accord with the change in the branch of the Liouville 
dressing. 



1. Introduction 



Large N matrix models [1] have proven to be a remarkable source of information 
about two-dimensional quantum gravity coupled to conformal matter with c < 1. They 
are the only available method for calculating sums over geometries to all orders in 
the genus expansion [2]. Some matrix model results have been reproduced directly 
in Liouville gravity [3-6], which gives us confidence that the discretized and continuum 
approaches describe the same theory. 

The matrix models which generate conventional discretized random surfaces are 
formulated with only single-trace terms in the action, such as Tr V(<&). Such models have 
been studied thoroughly, and their Liouville formulation is believed to be understood 
quite well. In gravitational dressing of any operator one encounters a two-fold ambiguity 
associated with the choice of branch of square root. For agreement with the conventional 
matrix models this ambiguity is always resolved by picking the branch which is smoothly 
connected with the semiclassical limit c — > — oo. Therefore, many results here are 
qualitatively semiclassical [7]. 

In this paper we study a different class of matrix models whose action, in addition to 
the single-trace terms, also contains trace-square terms such as g(Tr$ 4 ) 2 . Terms of this 
kind can glue a pair of random surfaces together at a plaquette. This contact may be 
thought of as a tiny neck (a wormhole), so that the network of touching surfaces may be 
assigned an overall genus and overall area. It is known that such microscopic wormholes 
are already abundant in the conventional theory with g = [8-10]. By changing g we 
are essentially changing the weight of some singular geometries in the path integral. It 
is not surprising, therefore, that a small increase in g does not change the universal 
properties of the model. If, however, we fine tune g to a finite positive value gt, then the 
universality class of the large area phase transition changes. This phenomenon was first 
observed [11] in a modified one-matrix model. It was found that there exists a critical 
value gt such that, for g < gt, the large area behavior of the sum over genus zero surfaces 
gives the string susceptibility exponent 7 = —1/2 characteristic of pure gravity, i.e. 

Fq(A) ~ A~ 3+ ^ ~ A' 1 ' 2 . 

For g > gt one finds branched polymer behavior with 7 = 1/2, which is not very 
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interesting because it corresponds to degenerate world sheets. Most interestingly, for 
g — gt there exists a new type of critical behavior with string susceptibility exponent 
1/3. This is the first example of a matrix model where new critical behavior occurs due 
to fine-tuned wormhole weights. 

Since ref. [11] a number of other such modified matrix models have been studied 
[12-16]. In general, as the trace-squared coupling is increased to a critical value gt, the 
string susceptibility exponent jumps from some negative value 7, found in a conventional 
matrix model, to a positive value 

7 = ^7- (1.1) 

7-1 

Essentially equivalent results have been obtained without using matrix models, on the 
basis of direct combinatorial analysis [17]. For a long time the positive values of string 
susceptibility exponent seemed very puzzling. Recently, however, a simple continuum 
explanation of these critical behaviors was proposed in ref. [16]. 

For all the conventional matrix models describing (p, q) minimal models coupled to 
gravity, the correct scaling follows from the Liouville interaction of the form 



A J d 2 aO m[n e a ^ , 



;i. 2 ) 



a+ = J—( v /l- c + 24h min - V25^~c) = - P+ L- 
2V3 V 2 W 

where O mm is the matter primary field of the lowest dimension, 



= 1 -&>-«)' (1 3) 

4pq 



A simple calculation reveals that the string susceptibility exponent is given by 



7 = 2 + ^-, (1.4) 



where Q = a/^V 2 - In ref. [16] it was argued that the effect of fine-tuning the touching 
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interaction is to replace the Liouville potential by 
A J d 2 aO min e a ^ , 



^( v / l-c + 24/ imin + v / 25^) = - P+ ll 1 



(1.5) 



Now the string susceptibility exponent is found to be 

7 = 2 + ^ = ^-. (1.6) 

a_ 7 — 1 

This establishes correspondence with the matrix model results, eq. (1.1). Thus, in the 
Liouville description of the modified matrix models we simply have to pick the other 
branch of square root in gravitational dressing compared to the conventional matrix 
models. This proposal has a number of interesting implications that are worth studying. 
For instance, using the scaling arguments of ref. [6] we find that the modified sum over 
surfaces of genus G should obey the scaling law 

d 2 F G 1 n 7) 

<9A 2 ~ A2G+7(1-G) ' { } 

In the matrix models this has only been checked for G = and 1. If true for all G, eq. 
(1.7) implies that it should be possible to define a double scaling limit for the modified 
matrix models. 

1.1. Summary of Results 

In this paper we carry out a non-perturbative study of various modified matrix 
models and confirm that eq. (1.7) is indeed valid. The simplest class of models we 
investigate are the modified multicritical one-matrix models [12,13], 



Tr V k ($) + (c 2 - A)Tr $ 4 - -L (Tr $ 4 ) ' 



;i.8) 



The critical potential of the fc-th model with g = is 

k 

vm = X>iy t+1 Q$ 2 * , 

1=1 

where q have been determined in ref. [2]. We choose to study the dependence of the 
sum over surfaces on a deformation of the potential by a term ~ $ 4 . For g — 0, it is 
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known that 7 = — \jk. The universal (leading singular) part of the sum over connected 
surfaces, F = logZ, is a function of the scaling variable t ~ (02 — X)N 2 ^ 2 ^\ In fact, 
for all g < gt the sum over surfaces is in the same universality class. For the fine-tuned 
value g = gt the string susceptibility exponent jumps to 7 = l/(k + 1). At this point 
there exists a critical value A c such that the universal part of the sum over connected 
surfaces, F, is a function of the scaling variable t ~ (A c — \)N 2 /( 2 ~i\ Our calculations 
reveal a remarkably simple non-perturbative relation 

00 

F(*) = log J dte ti+F ^ (1.9) 
— 00 

which connects the double scaling limit of modified matrix models with that of the 
conventional ones. We will show that this relation is quite general; it applies to all 
modified one-matrix and two-matrix models with trace-squared terms of the simplest 
kind, which describe c < 1 models coupled to gravity* 

There is a way, however, to alter the relation (1.9) by changing the type of trace- 
squared terms added to the action. We may replace the term which is the square of the 
lowest dimension operator by the square of some other scaling operator O. Adding gO 2 
to the action and fine-tuning g, we obtain a model where the gravitational dimension of 
the operator O changes from d to 

J=7-d (1.10) 

(the string susceptibility exponent remains unchanged). Remarkably, this change of 
dimension is again reproduced in Liouville theory by a mere change of the branch of 
gravitational dressing. Namely, if the operator is dressed by then 

a a 

where e a< ^ is the dressing of the Liouville potential. The relation (1.10) follows from 
(3+ + (3- = —Q. This gives us some further evidence in favor of the interpretation in ref. 
[16]. 



★ For c = 1 the relation is somewhat different due to the logarithmic scaling violations. 
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In the matrix models we, in fact, find the non-perturbative dependence on the cou- 
pling constant corresponding to O. If we perturb the action of the model by a term roO 
then, for g = 0, the universal free energy F is a function of t and to = toN 2 ^ 1 ^^^ 2 ^^. 
For g = g t we obtain a new universal free energy 

oo 

F(Mo)= log j dt e toio+F{t ' to) (1.11) 

— oo 

where to ~ ro7V 2 ( 1_d )/( 2_7 ). Thus, the calculation reduces to an integral over a different 
coupling constant from that in eq. (1.9). Now it is clear that, by simultaneously adding n 
different types of trace-squared terms and fine tuning their coefficients we can change the 
gravitational dimensions of n operators in the theory. The resulting partition function 
is a transform of the original partition function with respect to corresponding coupling 
constants U. The most general relation is 

F({*},{T})=logf] / dt ie ^=i tjtj+F{{tUT}) (1.12) 

where {T} is some set of other coupling constant that remain unintegrated. 

The organization of the rest of the paper is as follows. In section 2 we derive the 
integration over the coupling constants in eq. (1.9) via a trick familiar from wormhole 
physics. In section 3 we demonstrate the change of operator gravitational dimensions 
and derive eqs. (1.11) and (1.12). In section 4 we study the torus free energy which 
is directly related to the operator content of the theory. We attempt to interpret the 
modified matrix model results in terms of Liouville theory. We conclude in section 5. 
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2. Wormholes and Integration over Coupling Constants 



In this section we derive non-perturbative results for the sum over surfaces in mod- 
ified matrix models. 



2.1. The one-matrix models 



Our method is completely general but we will first illustrate it with the simplest 
example, the modified k = 2 one-matrix model which describes pure gravity. We need 
to study the following integral over an N x N hermitian matrix $, 



-TV 



Tr(i$ 2 -A$ 4 )-^(Tr<i> 4 ) S 



(2.1) 



A very helpful trick is to rewrite this as 



Z = 



N 



oo 

J dye~^~ J 



D$ e -^Tr(i$MA +J y)<I> 4 ) 



(2.2) 



This trick is generally useful when the action is perturbed by a square of an operator and 
has been applied extensively in wormhole physics [18]. In fact, the wormhole combina- 
torics leading to eq. (2.2) is the same as in four dimensions. Matrix models implement 
this combinatorics automatically, so that eq. (2.2) can be derived in one line. As in four 
dimensions we find integration over a coupling constant, which in this case is the quartic 
coupling. Here we are on a much firmer ground, however, because we know a great deal 
about the euclidean path integral, which is given by the logarithm of the matrix integral. 
This will allow us to obtain perfectly explicit and interesting results. 

The advantage of representation (2.2) is that we know how to perform the integral 
over $. The remaining integral over y is one-dimensional and can be thought of as an 
effective theory of baby universes, which lends itself to perturbative expansion around 
the saddle point. Relying on the well-known solution of the one-matrix model [19, 2], 



7 



we have 



log J V<S>e- NT *(^ 2 -( x+ y^) =N 2 (-a 1 x + ±a 2 x 2 ) + F(x,N 2 ) , 

F(x, N 2 ) = N 2 {-\a^ 2 + ...) + N°(-±- \ogx + . . .) + N~ 2 {a A x^l 2 + ...) + 0(N- 4 ) , 
x = c 2 - (A + y) . 

(2.3) 

We have chosen to separate the free energy into its singular part, F, and the leading 
non-singular parts of order iV 2 which play an important role in our discussion. From 
the leading order solution of ref. [19] we know the coefficients 

0,1 = 4 , ci2 = 576 , 03 = 6144^ , c 2 = 

4o . 

In the double scaling limit, t = xiV 4 / 5 ^ 5 is held fixed so that the subleading parts of 
F at each order in N become negligible. Thus, in this limit 

Fix, N 2 ) = Fit) = --t 5 / 2 - — log* + — r 5 / 2 + 0(r 10 ) (2.4) 

V , ) V J 5 24 6 216Q V ) V ) 



where we neglected x-independent additive terms. The expansion of F(t) follows from 

HP 



the fact that x(t) — ^dr satisfies the Painleve equation [2] 



d 2 X Q / + 2^ 



From eqs. (2.2) and (2.3), we have 



oc 




This may be thought of as an effective theory of baby universes. The interaction vertex 
of n baby universes arises from surfaces with n punctures (each puncture is generated 
by insertion of operator Tr<3> 4 ). The explicit quadratic term in x gives a kind of mass 
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term (inverse propagator) for the baby universe. Thus, the mass-squared is given by 
m? = i — a 2 , with the first contribution coming from the wormhole term in the matrix 
model, and the second from the degenerate sphere consisting of two plaquettes. We will 
analyze the three cases where the mass-squared is positive, negative and zero separately. 
The essential observation is that, since f(x) is of order N 2 , we may develop a large N 
expansion by integrating around the saddle point x s given by f'{x s ) = 0. 

Let us now show that, in the massive case g < l/a 2 the sum over surfaces is in the 
same universality class as the g = theory. Defining 



A = A c - A, 

we find that the location of the saddle point is 

A 



A c = c 2 - aig 



Xc 



1 - a 2 g 



+ C(A 3 / 2 ) 



Here the best way to analyze eq. (2.5) is by shifting the integration variable, z 
x — j_^ 2g . Discarding some non-singular terms in A, we have 



logZ(A,iV 2 ) = log / Ndzexp 



oo 

/ 



N 2 (l-a 2 g) 



z 2 + F z + 



N' 



2g \ 1 - a 2 g 

— oo 

After rescaling the variables, t = j^^N^^a 2 / 5 , z = zN 4 ^ 5 a^ 5 , we arrive at 

oo 

logZ(A,iV 2 ) = log J N^dzexp 



AT2/5 -4/ 5(1 _ } ■ 
-^z 2 + F{t + z) 



2g 



In the double-scaling limit, the gaussian term in the integrand becomes a delta-function. 
Therefore, the singular part of the sum over surfaces satisfies 

logZ(A,iV 2 )=F(t) , 

where F(t) is given in eq. (2.4). This explicitly establishes the universality for g < \ ja 2 . 
Thus, in their massive phase, the baby universes are irrelevant. 



★ In section 2.2 we show that this is true is general. The massive baby universes are irrelevant due 
to the general relation d > 7/2, where d is the gravitational dimension of the operator inserted by 
the baby universe, and 7 is the string susceptibility. In the case just studied, d — and 7 = —1/2. 
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For the tachyonic case g > 1/02, the saddle point near x = becomes unstable, but 
a stable saddle point appears at 



= £ + aA 1 / 2 + 0(A) , a>0, 



where x > is determined by the equation f"(x) = 0. We also deduce that f'(x) = 
N 2 A/g. Expanding f(x s ) in powers of A, we find that the leading singularity in the 
planar limit is 

logZ(A,iV 2 ) ~ iV 2 A 3/2 

which is indicative of the branched polymer phase. This behavior is quite generic because 
it does not depend on the precise form of F(x, N 2 ). 

The massless case g — l/a2 is critical. Here the position of the saddle point acquires 
a new scaling, 

Integrating around the saddle point in eq. (2.5), we find 



hgZ 



+ C(l/iV 4 ). (2.6) 



We find that all the terms in this expansion are important. After some calculation, we 
arrive at 

log Z = N 2 (- A 5 / 3 + ...) + iV°(-— log A + . . .) + N- 2 ( —A" 5 / 3 + ...) + 0(N- 4 ) , 
& v 5 ; v 36 ; v 960 ; v ; 

A = Aa2G 3 2//5 . 

Thus, the singularity for any genus occurs as A — > 0. The structure of the leading 
singular terms suggests that we may now define double scaling limit by keeping the 
variable t ~ AN 6 / 5 fixed. In fact, this follows directly from eq. (2.5). For g = l/a,2 the 
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singular part of the sum over surfaces, given by F = log Z, becomes 

oo 

F(A,iV 2 ) = log J dxNe N ^ Ax+F ^ N ^ . (2.7) 

— oo 

Introducing scaling variables 

t = xN^af , i=AN 6 / 5 a 2 a 3 2/5 , 

we find that 

oo 

F(A, N 2 ) = F(t) = log J dte a+F ® . (2.8) 

— oo 

This is our main result, which establishes a simple relation between the double-scaling 
limits in the modified and conventional matrix models. It is uncertain whether eq. (2.8) 
has a truly non-perturbative meaning: there are well-known problems in defining F(t) 
non-perturbatively. Even if they are overcome, it is not clear if the integral over t will 
converge. What is certain is that eq. (2.8) determines the sum over modified surfaces 
of any genus, i.e. it works to all orders of perturbation theory. Using saddle-point 
techniques, summarized in eqs. (2.5) and (2.6), we may generate the large t expansion 
of F directly from the integral representation (2.8) and the large t expansion of F. 

Eq. (2.8) applies equally well to matrix models with non-symmetric potentials, which 
are more basic because they do not double the degrees of freedom. For such models we 
divide F(t) from eq. (2.4) by 2 to obtain the sum over surfaces. After a redefinition of 
t, we find 

F(t) = _| t v 2 __L logt + _l_ t -v 2 + 0(t - I o ) . 

Substituting this into (2.8), and generating the saddle-point expansion, we find the 
modified sum over surfaces for non-symmetric models, 

F(t) = \t>'* - I log* + |^f-*/3 + O (f-10/3) 

where we have explicitly calculated the contributions up to genus 2. 
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Extension of the methods presented above to the k-ih multicritical one-matrix model 
is quite straightforward. First we rewrite eq. (1.8) as 



Zu = 



— oo 

Using the variable x = c 2 — A — y, we find 



J dye'^ Jv<S>e~ NT < v ^+^- x -y^ 4 ) (2.9) 



oo 

N 



dxe fk{x) , 



-oo 



N 2 r o 9 n 9 (2-10) 

fk(x) = — [(ca - A) + 2x(c 2 - a l5 - A) - x (1 - a 2 g)] + F fc (x, N 2 ) , 
F fe (x,iV 2 ) = -^ I a 3 ivV 2fc+1 )/ fc + ... 

In the double scaling limit of the g — theory, t ~ x ]\f 2k /( 2k + 1 ) is held fixed, and the 
sum over surfaces is given by F^{t). Analysis of the saddle-point expansion shows that 
the theory is in the same universality class for any g < l/a 2 . 

For g — l/a 2 , we instead have 



oo 

F fc (A,iV 2 )=log J dxNe N2a ^ x+F *^ N2 \ 



A = c 2 A . 

a 2 

Introducing scaling variables 

t = xN 2k /( 2k +^4 /{2k+1) , t = AM 2k+2 ^ 2k +^a 2 a~ k/{2k+1) 

we arrive at the modified sum over surfaces in the double scaling limit, 

oo 

F k (t) = log J dte tl+Fk ^. (2.11) 

— oo 

From the fact that the modified sum over surfaces is a function of t, and that the 
original expansion was in powers of l/N 2 , it follows that the genus G contribution scales 
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(2fc + l)(l-G) 

as t k + 1 . Using the fact that, for non-symmetric matrix potentials, 

F k (t) = -^-/ 2k+1)/k -^lo g t + j^ aj (k)t-^» k (2.12) 

i=i 

we generate the genus expansion of with the saddle-point methods, 

— .—. k+l^ i 1 fk-1 l\ - o, 1 N 

This confirms the known result [12, 13] that on a sphere the string susceptibility exponent 
is 7 = fcqrp Thus, the order of the phase transition for planar surfaces has changed from 
the third to the second. We have also established that the susceptibility exponent at 
genus G is 7 + G(2 — 7), in agreement with eq. (1.7). Since this genus dependence has 
such a natural explanation in Liouville theory, it provides a solid argument in favor of 
the Liouville interpretation of the modified matrix models [16]. 

2.2. The two-matrix models 

In this section we consider a more general class of matrix models and show that, with 
the simplest kind of trace-squared terms, the non-perturbative relation (1.9) applies to 
them as well. 

First we address the Ising model coupled to gravity, which is well known to be 
described by a two- matrix model [20]. Its simplest modified version was introduced in 
ref. [16], 

Rising = / V^ l V^ 2 e- N ^+ S ^+ kT ^^ , 

J (2.13) 
S{$) = Tr ($ 2 - A$ 4 ) - -| (Tr$ 4 ) 2 . 

We have added trace-squared terms of the form 

s[(Tr$f) 2 +(Tr$ 4 ) 2 ] (2.14) 
which implies that the value of the Ising spin at the two ends of a wormhole is required 
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to be the same* Rewriting eq. (2.14) as 



f (Tr($? + $t)) 2 + f (Tr($t-$|)) 2 (2.15) 
and applying our trick to each of the two terms, we arrive at 

oo 

Rising = 7T- / d y dve 29 / V<S>!V<S> 2 e~ NS , / n 

2*9 J J (2.16) 

— 00 

S = Tr [$? + $l + A;Tr$i$2-(A + y)($i + $|)-t;($f-$|)] . 

We perform the matrix integral first, and save the integrals over v and y until the end. 
The matrix integral describes the Ising model in magnetic field v [20]. If we tune k to 
its critical value and define x = c 2 — A — y, where c 2 is the critical quartic coupling, then 

log J £>$i£>$ 2 e~ NS = N 2 {-aix + ^a 2 x 2 ) + F (x, N 2 ) + F^x, v, N 2 ) , 

F (x, N 2 ) = N 2 (~a 3 x 7 / 3 + ...) + N°(~ \ogx + ...) + 0(N~ 2 ) , (217) 

oo 2n 

71=1 

The v 2n vertex in the diagrammatic expansion for v is given by the connected correlation 
function of 2n operators Tr (<&f — $2). This operator is the gravitationally dressed spin 
field, which is known to have gravitational dimension 1/6. Therefore, 

(Tr ($1 - $|)) 2 ) = N 2 (b + b'x 2 / 3 + ...)+ iVV 5 / 3 + . . . 



The inverse propagator (mass-squared) for v is then given by Q — bj. By an explicit 
calculation, following the results of ref. [20], we determine 

84024 , 41256 

a? = , = . 

1 625 ' 625 

It is important that a 2 > b. Let us imagine dialing g up. For g < \jai both the x and 
v baby universes are massive and contribute only subleading terms to the free energy. 



★ In the next section we will relax this condition and allow a spin flip when a wormhole is traversed. 
We will see that more general theories can be constructed this way. 
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Here we find the same universality class as the unmodified Ising model. For g — 1/02, x 
becomes massless and changes the critical behavior, v is still massive, however, so that 
the integral over v has the form 



Ndv exp 



_IjVV(a 2 - 6) + F^ niveTS& \x, v, N 2 ) 



In terms of the rescaled variables v = vN 2 ^ 1 d ^>^ 2 ^ and t ~ xN 2 ^ 2 7 \ this becomes 



/ 



2d-7 

N 2 -T dv exp 



1 4d-2 7 „ 

--N v 2 (a 2 -b)+ Fi(t,v) 



If 2d — 7 > 0, then in the double scaling limit the gaussian factor becomes a delta 
function, so that v is frozen at zero* Since Fi(t, 0) = 0, the v integral does not contribute 
to the effective action for x. Thus, all the dominant terms in Z\ s ^ can be calculated 
from 



N 



sing 



dxe f W , 



(2.18) 



N 2 



f(x) = TT [( C2 - A ) 2 + 2x ( C2 ~ ai9 - A) - x 2 (l - a 2 g)] + F (x, N 2 ) 



From here on the calculation is analogous to those encountered in the one-matrix model. 
Introducing scaling variables 



t~(c 2 -^-A)iV 8 / 7 , 

02 



t ~ xN 6 / 7 



we arrive at the relation (1.9). 



★ This holds here because the gravitational dimension of the spin operator is d = 1/6 and 7 = —1/3. 
In fact, this holds in general because in all conventional matrix models d > and 7 < 0. We 
conclude that all massive baby universes are irrelevant. 
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The calculation presented above for the Ising model can be carried over almost 
verbatim to any modified two-matrix model of the form 



Z= f V$ 1 V$ 2 e- N[Sp (* l)+s ''(* 2 ' )+ '* rr * 1 * 2] , 

J ' ' (2.19) 

5 p ($) = Tr ($ 2 - A$ 4 + . . . + r$ 2p - 2 ) - -| (Tr $ 4 ) 2 . 



The parameters of the potential S p (&) can be tuned [26] in such a way that this describes 
an arbitrary minimal model coupled to gravity. For any such model, the modified sum 
over surfaces F(t) is related to the conventional sum F(t) by eq. (1.9). The general 
validity of this relation implies the generality of eq. (1.1). To show this, note that 



Since tt ~ AxN 2 , we have 



1 + ^ = i 



2-7 2-7 

from which eq. (1.1) follows. Thus, if the asymptotic expansion of F(t) is in powers of 
t 2 ~ 7 , then the asymptotic expansion of F(t) is in powers of . 



2.3. c 



One interesting theory that remains to be discussed is the c = 1 model coupled to 
gravity. We will consider compact target space of radius R, which is described by matrix 
quantum mechanics at finite temperature [21]. The path integral that generates the sum 
over touching surfaces is [14,15] 



= / P$(t) e Jo 



Tr ( i $ 2 + \ $ 2 - A$ 3 ) - ^ (Tr $ 3 ) ' 
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with $(27ri?) = $(0). Let us introduce the normal mode operators 

2ttR 

P= J dtTr$ 3 (t) , 



2nR 

C n = -j= J dtcos-Tr$ 3 (t) , 



2nR 

If nt 
Sn = ^J dt sin— Tr$ 3 (t) , 



and write the trace-squared term as a sum of squares 

2 « R / oo 

The operators C n and S n are known to have gravitational dimension d = n/2R [22]. 
Using the by now familiar trick, we introduce a "baby universe variable" for each squared 
operator in the action to derive 

Z~ j dy H dy n dz n e'^ { y l+y - +zl) 

-oo n=1 

j V<$>{t) e- N C Rdt [Tr(i<& 2 + i<l> 2 -(A+ y o)$ 3 -E l C : i (y^ l+ , l S l ))] 

Performing the matrix integral, we get 

log J V$(t) e~ NS = 2TvRN 2 (- ai x + l -a 2 x 2 ) + F (x, N 2 ) + F^x, y n , z n , N 2 ) , 



F (x, N z ) = RN\-a^/ log* + ...)- —(R + -) log* + . . . (2 . 20 ) 

oo 

F^ymZn.N 2 ) = nRN 2 J^ivl + 4)% + &nO*/l logx\) n / R + ...) + ••• 



n=l 



We have exhibited the terms in F\ that come from the two-point functions of C n and 
S n . These quadratic terms determine whether the variables y n and z n become critical 
simultaneously with the variable x = c 2 — A — yo. 
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From a calculation of the momentum dependence of the puncture two-point function 
[22], we have 



oo 

„2 



-oo 



where f(s) is proportional to the Fourier transform of the classical trajectory at the top 
of the critical potential. Since b n is a decreasing function of n, and &o = °2, we conclude 
that b n < CL2 for all n > 1. This crucial finding implies that, as the variable x becomes 
critical for g = 1/02, all the other baby universe variables are still away from criticality. 
Since their propagators are massive, their fate is the same as of the variable v in the 
Ising case: integrating them out makes no effect on the relevant terms in the effective 
action for x. For g = 1/02, the important integral over x reduces to 

00 

F(A,iV 2 )=log J dxNe ^RN 2 Ax+F ( X ,N*) (2 21) 

—00 

where A = 02(02 — — A). While in other models we could express this integral directly 
in terms of scaling variables, for c = 1 this is impossible because of the logarithmic 
scaling violations in Fq(x,N 2 ). Actually, the situation turns out to be even simpler 
than for c < 1. 

In the leading saddle point approximation, 

F(A, N 2 ) = 2nRN 2 Ax s + F (x s , N 2 ) , 

_ m = = 2M (2.22) 

ox 

Thus, F is the Legendre transform of —Fq, with 2ttRN 2 A being the conjugate variable 
of x. The leading order relation between A and x is 

A log A ~ x , 

in agreement with ref. [15]. Our analysis of integration around the saddle point, based 
on eq. (2.6), indicates that, remarkably, all such corrections are suppressed by powers of 
log A. Thus, in the double scaling limit where iVA is kept fixed, eq. (2.22) is exact. This 
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Legendre transform was introduced in ref. [21] to calculate the sum over "one puncture 
irreducible" surfaces. It was shown to satisfy a simple equation, 



d 2 F 

M2 



2nRN 2 p(A) = RN" 1 



lnA+]T (2NAVR) 2m f m (R) 



m=l 



/ m ( J R)-(2m-l)!^|2 -2||2 " 2 I (2fc)![2(m _ » 



(2.23) 



fc=0 



where p(A) is the temperature corrected density of states, and A is the distance of the 
Fermi level from the top of the potential. Integrating eq. (2.23), we find 



il-(2iVAv / ^) 2 lnA-2/i(i?) lnA+ V fm+1 ^ (2NAVRY 
8 I ' v ' ^ m(2m + l) v ' 

\ m=l v ' 



■2m 



, (2.24) 



where we have exhibited only the terms that survive in the double-scaling limit. It is 
remarkable that the c = 1 model with fine-tuned wormhole weights directly generates 
the "one puncture irreducible" surfaces. This model, which has no scaling violations as a 
function of the area [15], is in many ways simpler and more natural than the conventional 
c = 1 model. 



3. New Gravitational Dimensions 

One lesson we can draw from the preceding section is that, even though a given 
model may have many baby universe integration variables, it is usually the case that 
only the integral over the lowest dimension coupling affects the critical behavior. It is 
clear, however, that this cannot be the most general situation. In this section we show 
how to make other integrations relevant by changing the type of trace-squared terms 
added to the action. 

As an instructive example, let us consider the modified Ising model of section 2.2 
with a more general class of trace-squared terms, 

|(Tr($t + $|)) 2 + |(Tr($t-$|)) 2 . (3.1) 

For g 7^ g' this introduces a term of the form Tr <3>^Tr $| which generates wormholes 
with opposite values of spin at the two ends. It is not surprising that such processes 
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can make the integration over the spin field coupling constant relevant. We may, for 
instance, set g — (actually, any g < l/a2 will do), while fine tuning g' to its critical 
value. The partition function becomes 



J J ' (3.2) 

— oo 

S = Tr [$? + $| + A;Tr $!$ 2 - A($f + $|) - (v + r a )($f - $|)] 

where we have introduced coupling constant r a in order to study correlation functions 
of Tr ($| — $|). Defining a shifted variable u = v + r a , we perform the matrix integral 
first and reduce the modified free energy to 



F = F (t) +log J due f ^ , 



N 2 

/(«) = "^(« 2 - 2«r a + r a 2 ) + F X (A, u, iV 2 ) , (3-3) 
Fi = -N 2 u\b + b'A 2d -~< + ...) + 7V 2 w 4 6"A 4d - 2 -T + . . . 

where t ~ AN 2 ^ 2 ^ and A = C2 — A. The universal part of F\ is a function of t and 
t a = mA^ 2 *^^^^ 2-7 ). In our specific case 7 = —1/3 and d = 1/6. If we now fine tune 
g' — 1/b and introduce the scaling variable 



t a = br ff N 2 ^-^ 2 -^ , 
d — 7 — d , 



(3.4) 



then the universal part of the modified free energy is given by 

00 

F(t,i a )= log J dt a e t ^ +F ^ (3.5) 

—00 

Here F(t,t a ) is the universal part of the conventional sum over surfaces. 

Eq. (3.4) implies that the gravitational dimension of the spin field has changed 
from d = 1/6 to d = —1/2. Thus, in modified matrix models negative dimensions arise 
naturally. Although we have discussed a specific example, the change of gravitational 
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dimension given by eq. (3.4) is general. As shown in section 1.1, this formula agrees with 
the change in dimension caused by changing the branch of Liouville dressing. Therefore, 
there are serious reasons to believe that such operators, which were previously thought 
not to exist, are in fact present in the spectra of modified matrix models. 

Eq. (3.5) shows that, by a fine-tuning of g', the coupling constant corresponding to 
the spin field has been driven to criticality. It is not hard to see that a simultaneous 
tuning of g to l/a2 also makes the coupling constant t, corresponding to the puncture 
operator, critical, so that 

oo 

F(i,t a )= log J dtdt a e ti+tJ ° +F ^ . (3.6) 

oo 

It is now clear that a fine-tuning of n parameters in the trace-squared terms can result 
in integration over n coupling constants, giving the general formula (1.12). 



4. Sum over Surfaces of Genus One 

In this section we focus on the torus contribution to the free energy which, in any 
string theory, is directly related to the spectrum. For all conventional matrix models 
the torus free energy has been successfully reproduced by path integration in Liouville 
theory [23,24]. For this reason it is particularly interesting to calculate the corresponding 
quantity in modified matrix models and ask for its continuum interpretation. 

For any (p, q) minimal model coupled to gravity the torus free energy is* 

This result was reproduced [23] in Liouville theory with the interaction term of eq. (1.2). 
Let us now study the modifications to this result due to the integration over coupling 
constants. Consider, for instance, a model where the gravitational dimension of operator 
O has been changed from do to 7 — do- Here the modified sum over surfaces is given 



★ We are quoting the answer for matrix models with non-symmetric potentials. For corresponding 
models with symmetric potentials the free energy is doubled. 
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by eq. (1.11). Setting to = and performing gaussian integration around the saddle 
point, we find 

F G =\t) = F G =\t) + i( 7 - 2d Q ) log* . (4.2) 

This result is puzzling from the point of view of the simplest Liouville approach. Since 
we have not changed the string susceptibility exponent, it would seem that the Liouville 
action is still (1.2), and that the calculation of ref. [23] with the result (4.1) should still 
apply. The matrix model tells us otherwise: the moment we change the gravitational 
dimension of an operator, the torus free energy receives a correction. We may speculate 
that in Liouville theory this correction originates from a boundary term in the modular 
integral, but at the moment we do not know how to derive it directly. In the following, 
however, we will give a plausibility argument for the presence of the correction found in 
eq. (4.2). 

Our argument is based on the interpretation [25] of the torus free energy in (p, q) 
models as the sum over zero-point energies of an infinite number of one-dimensional 
particles (harmonic oscillators). Each oscillator corresponds to an operator in Liouville 
theory of the form Oe^, where O has dimension h and 



co = ^l^-2 + 2h. 

uj, which is the "Liouville energy" , gives the frequency of the oscillator. A priori there is 
a sign ambiguity for u, but in conventional Liouville theory all ui are taken to be positive. 
Each operator contributes zero-point energy ^oo to the coefficient of the Liouville volume, 
— \ogt/a, where a = a + from eq. (1.2). Thus, 

^'w = -^Ef («) 



a 



For the (p, q) model the spectrum of energies is given by [25] 



Si 



y/2pq 

where Si are the positive integers not divisible by either p or q. After substituting this 
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into (4.3) and using zeta function regularization of the infinite sum, ref. [25] recovered 
eq. (4.1). 

When we fine tune the theory as in section 4, we replace uo by —too in the Liouville 
dressing of operator O. If we make the same replacement in eq. (4.3), we arrive at 

F G=1 (t) = F G=1 (t) + ^ log t . (4.4) 

a 

Remarkably, since 7 — 2do = this agrees with the matrix model result, eq. (4.2)! 

This suggests a connection between the operator content and the torus free energy in 
modified matrix models. Our basic premise is that a fine-tuned operator with negative 
Liouville energy, — contributes a negative zero-point energy, — \\to\. This suggests 
that, from the space-time point of view, a negatively dressed operator is a fermion. It 
would be interesting to find an explanation for this effect. 

Proceeding to other modified matrix models we note that, for the models described 
by relation (1.9), 

F G=1 (t) = - fr-^-^ logt ^ log*. (4.5) 

1 ; 24(p + g + l) S p + q + 1 S K J 

A naive Liouville calculation with potential (1.5) would give only the first term in the 
above [16]. However, applying eq. (4.3) with a = a_ from eq. (1.5), and including 
negative zero-point energy for O m [ n , reproduces the matrix model result. 

The second terms in eqs. (4.2) and (4.5) are due to integration around the saddle 
point in eqs. (1.11) and (1.9). These corrections can be eliminated if we treat the baby 
universe variables classically, i.e. if we freeze them at their saddle point values. In such 
a theory F would simply be the Legendre transform of —F. Of course, this theory does 
not correspond to the original matrix model with trace-squared terms, but it does have 
a simple geometric interpretation. It calculates the sum over trees of touching random 
surfaces (bubbles), with each bubble allowed to have arbitrary genus. In other words, 
the wormholes are present, but they are not allowed to increase the overall genus. Since 
from the world sheet point of view this constraint is highly non-local, we do not regard 
such a theory as natural. Eq. (4.5) shows, however, that the relative importance of 
surfaces where a wormhole closes the loop decreases with increasing p and q. For c = 1 
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the quantum effects associated with the baby universe variables become completely 
negligible, so that the modified sum over surfaces is simply the Legendre transform, eq. 
(2.22). 

5. Conclusions 

Our non-perturbative solutions of matrix models modified by various trace-squared 
terms strongly suggest that there exists a continuum Liouville formulation of these mod- 
els. All the modified scaling exponents agree with the idea [16] that a fine-tuning of 
trace-squared terms changes the branch of Liouville dressing of some operators. Since 
the new branch of dressing does not have a semiclassical limit, the resulting Liouville 
theory is more complicated and more interesting than the conventional one. We hope 
that our matrix model results will provide a useful guide towards such a theory. 

The solution of the modified matrix models is also quite interesting in itself because of 
its connection with general wormhole phenomena in quantum gravity. The microscopic 
wormholes, introduced by the trace-squared terms, lead to integration over coupling 
constants, as expected on general grounds. Such integration arises in any theory with 
bilocal operators in the action. Physical effects of integration over coupling constants 
have even been found in theory of elasticity, where they change the order of phase 
transitions* [27]. Our work provides another example of a system where integration 
over coupling constants introduces a profound change, affecting even the order of the 
phase transition for planar surfaces. It is interesting to look for other physical systems 
where coupling constants turn into dynamical variables. 
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